Development of accurate solutions for a 
classical oscillator 



Paolo Amore 

Facultad de Ciencias, Universidad de Colima, 
Bernal Diaz del Castillo 340, Colima, Colima, 
Mexico. 

Nestor E. Sanchez 

Department of Mechanical Engineering, The University of Texas at San Antonio, 

San Antonio, TX 78249 



Abstract 

We present a method to obtain arbitrarily accurate solutions for conservative clas- 
sical oscillators. The method that we propose here works both for small and large 
nonlinearities and provides simple analytical approximations. A comparison with 
the standard Lindstedt-Poincare method is presented, from which the advantages 
of our method are clear. 

Key words: Nonlinear oscillator, Lindstedt-Poincare 



1 Introduction 



In this paper we consider the problem of calculating the periodic solutions to 
the differential equation 

m^ + /(x)=0. (1) 



Eq.(l) is the Newton equation for a particle of mass m moving under the ac- 
tion of a force f{x). We assume that the system is conservative and therefore 
we introduce the potential V{x) = — J f{x)dx. The exact solution of Eq.(l) 
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is possible only in a few cases: more often one only disposes of approximate 
results, which usually arc obtained by applying pcrturbative methods. Per- 
turbation methods arc based on an expansion in some small parameter in the 
problem and the approximate solutions are therefore obtained in the form of 
a polynomial in such a parameter. This is the case of the Lindstedt-Poincare 
method and of the multiple-scale method, which are widely used in the com- 
munity. 

Unfortunately the validity of the perturbative approaches is restricted to the 
domain of small parameters and the series obtained in such a manner have a 
finite radius of convergence: in other words these series become useless when 
the parameters are larger than the radius of convergence. While other tech- 
niques are been developed in the literature to deal with this problem, see for 
example [1,2], we wish to present a novel method which was recently devised 
by one of us [6,7,8] and which allows one to calculate accurate analytical so- 
lutions for a classical oscillator, described by Eq. (1). The method is based 
on the powerful ideas of the Linear Delta Expansion (LDE) [3], which allows 
to obtain, to any given order, fully analytical and extremely accurate expres- 
sions. Since our method is not based on an expansion in a small parameter, 
the series obtained converges regardless of the values taken by the parameter 



Since our method is both systematic and analytic it provides several advan- 
tages with respect to other techniques which also claim to work for arbitrary 
functions f{x): this is the case for example of the method of He [4], whose 
shortcomings have been evidentiated by one of us in [5]. 



2 The method 

This problem has been previously studied by one of us in [6,7], where a non- 
perturbative method which allows to calculate the period of the oscillations 
was devised. 

We will now briefly review the method and then show how to generalize it 
to calculate the solutions. For a conservative system the total energy E is 
constant and the period is simply given by 



where x± are the inversion points, obtained by solving the equation E ~ 
V{x±). Notice that V{x) = — J f{x)dx is the potential: 



itself. 




(2) 
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As discussed in [6] the problem of calculating this integral can be converted 
to 

Ts= [ '^'^ dx, (3) 

Eo - Vo{x) + S[E -Eo- V{x) + Vo{x)] 



where Vo{x) is a potential chosen to interpolate the original one. It is expected 
that Vo{x) is simpler than V{x). For S — 1 Eq. (3) reduces to Eq. (2), whereas 
ior 5 — this formula yields the period of oscillation in the potential Vo{x). We 
will treat the term proportional to 5 as a perturbation and expand in powers 
of 6. Since Vo{x) depends upon one or more arbitrary parameters (which we 
will indicate with A) a residual dependence upon these parameters shows up 
in the period when the expansion is carried out to a finite order. In order 
to eliminate such unnatural dependence we impose the Principle of Minimal 
Sensitivity (PMS) [10] by requiring that 

dT 

ax = °- 



We now proceed to define 



_ E-Eq- V{x) + Vojx) 
^^""^ ~ Eo - Vo{x) ■ 

and write Ts as 

n=0 ^- ^ ^^Eo - Vo{x) 



provided that the series converges uniformly, which is the case if |A(a;)| < 1 
for every x, x^ < x < x^. We can choose Vo{x) so that each integral in (6) 
can be performed, thus obtaining an equivalent series representation for the 
original integral (2). 

Amore and collaborators have applied this method to obtain very precise 
analytical approximations for the period of several oscillator, showing that 
the error decreases exponentially with the order of the approximation ^ . 

In particular they considered the Duffing oscillator, which corresponds to the 
potential V{x) = ^ x'^ + j x"^. The interpolating potential was chosen to be 

^ Recently, Amore and Arceo have given in [8] a precise estimate for the rate of 
convergence of the series (6) and they have shown that it converges in ah the physical 
region. 
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Vo{x) = -^x% where A is an arbitrary parameter. Working to first order 
they obtained the optimal value \pMS — \/^A/2 and the simple formula: 



which provides an error less than 2.2% to the exact period for any value of n 
and A. 

Physically the PMS allows us to obtain the best potential around which to 
perform our expansion: the form of the potential depends upon the arbitrary 
parameter A and will in general depend on the order to which the calculation 
is made. 

Wc will now generalize the results of [6,7] by using the conservation of the 
energy to obtain the solution to Eq.(l). It is straighforward to see that the 
time spent by the body to go from x- to a point X [x- < X < Xj^) is given 

by 

^ dx . (8) 



After repeating the procedure previously explained we have 



».o "! 2" I - V„(x) 

The optimal parameter A will again be chosen according to the same criteria 
previously adopted. Notice that, once the integrals in (9) are calculated, one 
obtains a convergent series representation for the time i as a function of the 
position X: having proved in [6] the convergence of (9) for X — x+, the 
convergence for X < x+ follows. 

As an application we consider the Duffing oscillator and working to first order 
we obtain 



^ ^ (6M^ + 8) cos-^ (f ) - A^X 

(3^2 + 4)3/2 

or equivalently 
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PMS 



t = COS 



-1 fl'^ _ ^'"^ ^ 



A2 



+ 8 



In Fig. 1 we have compared the numerical solution corresponding to fxA^ — 10^ 

with the approximation given by Eq. (11): otu^ simple analytical formula pro- 
vides an excellent approximation, even in presence of a huge nonlinearity. This 
regime is clearly outside the region of applicability of perturbative methods, 
such as the Lindstedt-Poincare (LP) method. 

To ease the comparison with the LP method we notice that the solution to 
the nonlinear equation will have the general form 



X{t) = 5Z Cn cos [{2n + l)nt] . 



n=0 



We can obtain the Fourier coefficients c„ simply by using Eq. (11): 



(12) 



2 f at 

c„ = / ^ TV cos[{2n + l)npMst] dX 

TT J aX 

-IjiA^ + 2/iX2 - 8 



= — X 

71 J 
-A 



COS < (2n + 1) 



2^(3^2 + 4^ V^-?. 

^ fX\ A^iX ^Jl-§ 



cos 



a) 



+ g 



dX 



(13) 



Eq. (13) cannot be evaluated analytically in its present form. However we 
notice that the function 



^^_Ai,X^ 



A2 



6/iA2 + 8 



(14) 



fuUfills the constraint |^| < —1/12 and therefore can be used as an expansion 
parameter. 

Working to order we obtain 



SPMS) _ A (2644 + 107456/12^4 + 145408/1^2 + 65536) 
'° ~ 1024 (3//A2 + 4)^ 

(PMS) _ (6435/13^6 + 26424/12^4 + 36096/1 A^ + 16384) 



Ci 



2048 (3/1^2 + 4)' 



(15) 
(16) 
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(PMS) 
'^2 



(PMS) 



5^V^ (427/i^>l^ + nU^A^ + 768) 

6144 (3M^ + 4)^ 
49AV^ (5/uA2 + 16) 
12288 (3/x A2 + 4)^ 



(17) 
(18) 



Notice that the coefficients above are rational functions of fiA^ and all tend 
to finite values for /lA^ — > 00. For /lA^ -C 1 the perturbative expressions are 
found 



JPMS) 



SPMS) 



(PMS) 
'-2 



^{PMS) 



A 

A 

A 
A 



IjA" 51/^2^^ ^ 1287At=^A6 ^ 



(M^)^ + . . . 



32 

5 

2048 
49 



2048 

{^Ar - 



65536 



745 
196608 



196608 



(M^)^ + ■ ■ ■ 



(19) 
(20) 

(21) 
(22) 



which in part reproduce the results obtained with the LP method: 



X^^^\t) = A 1 - ^ + — ^- ^Jl^^ + ... cos pLP t] 



+ A 



+ A 



+ A 



32 1024 32768 
3/12^4 29/X-M6 



. 32 
^1024 
32768 



128 

•i^i^A^ 



16384 
+ . . . cos [5Qip 



cos [3f2Lp t] 



2048 

+ . . . ) cos [IQlp t\ + ... 



(23) 



It is useful to consider the exact solution to the Duffing equation which, given 
the initial conditions used can be cast in the form 



X(--)(t)^Acn(yrT^t|^^^) . 



(24) 



We can use Eq. (16.23.2) of [9] to write 
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cn(u\m) = —;= — ; — r — — cos(2n + l)f (25) 

where q = e-^-^(i-'")/^M and v = 7ru/2K{m). 

We can easily read off Eq. (25) the Fourier coefficients of the cn function 
and compare them with the approximations Eqs.(15),(16),(17) and (18): the 
leading Fourier coefficient is reproduced with a maximum error of 0.15%. In 
Fig. 2 wc plot the absolut value of the error S = l-^ — ^exac" — ) x 100. We 
notice that for fiA^ < 1 the error is quite small and dies exponentially fast as 
the limit — > is approached; on the other hand, in the limit iJ,A^ — > oo the 
error reaches a plateau. The plateau exists because our approximate coefficient 
^approx j^^g ^j^g correct asymptotic behaviour, i.e. 



^ -77^ ^ 1-0017 . (26) 



M^^oo cr-^ 110592 r(|)^ 



Although the remaining coefficients are reproduced with less accuracy (for the 
coefficient Ci we have a maximum 10% of error) the overall solution is very 
accurate, since Cq is much larger than all the remaining coefficients. To prove 
this statement we can read off the exact expressions for the c^^"'^* from Eq. (25) 
and calculate 



„exact pUTT /I I p7r\ 

^ lim ^ = ^ , ^ ' , (27) 



which decays exponentially for large n. In the case n = 1, we see that, even in 
the asymptotic limit jiA^ oo, the coefficient Cq^"^* is about 22 times larger 
than c^^"'^*, which is the key to the precision of our results. 

Notice that the limit /lA'^ oo corresponds to considering a purely anhar- 
monic oscillator V{x) — | x^: clearly, our method is capable to deal quite 
efficiently also with this case. 

It is worth stressing that the present analysis has been carried out only to first 
order 5: however, since the method is geometrically convergent ( see [8]) one 
expects that much higher precision can be obtained by applying it to higher 
orders, although this issue is not pursued in this paper. We also stress that 
our method is completely general and that it can be applyied to a large class 
of potentials: a detailed analysis of the application of the method to calculate 
the period of general oscillators is given in [7] . 
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3 Conclusions 



We have presented a method to obtain arbitrarily accurate solution for a con- 
servative oscillator. The particular technique that we used worked for small 
and large nonlinearities of the equations, and provides a simple but very ac- 
curate approximation. A comparison with the exact solution and with the 
perturbative solution obtained using the Lindstedt-Poincare is done. It was 
shown that errors as small as 0.15% are recorded on the leading Fourier coef- 
ficient. 

The method provides several advantages over other well established methods 
in the literature: first of all, to the best of our knowledge our method is the only 
nonperturbative method which allows to obtain fully analytical results and for 
which the exponential convergence of the series is proved ([6,7,8]); secondly, 
previous work done in [6,7,8], where the method was used to calculate the 
period (and not the solution) has shown that our method can be used for a 
quite large of class of potentials, even in cases where the exact result cannot 
be obtained; finally it is easy to calculate higher order contributions {never 
involving special functions) with our method. 
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Fig. 1. Solution corresponding to Eq. (10) for /lA^ = 10^ (solid curve). The dashed 
curve is numerical. 



10 



0.001 



[i] le-06 - 



le-09 



le-12 
0.0001 




10000 



Fig. 2. Absolute value of the error over the Fourier coefficient cq using the approxi- 
mate expression (15). 
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